US008995724B2 


d2) United States Patent 

Asner 


(io) Patent No.: 

(45) Date of Patent: 


US 8,995,724 B2 

Mar. 31, 2015 


(54) REMOTE SENSING ANALYSIS OF FOREST 
DISTURBANCES 

(71) Applicant: Carnegie Institution of Washington, 

Washington, DC (US) 

(72) Inventor: Gregory P. Asner, Half Moon Bay, CA 

(US) 

(73) Assignee: The Carnegie Institution of 

Washington, Washington, DC (US) 

( * ) Notice: Subject to any disclaimer, the term of this 

patent is extended or adjusted under 35 
U.S.C. 154(b) by 110 days. 

(21) Appl.No.: 13/800,726 

(22) Filed: Mar. 13, 2013 

(65) Prior Publication Data 

US 2013/0216103 A1 Aug. 22, 2013 

Related U.S. Application Data 

(63) Continuation of application No. 13/459,804, filed on 
Apr. 30, 2012, now Pat. No. 8,416,994, which is a 
continuation of application No. 12/083,833, filed as 
application No. PCT/US2006/040814 on Oct. 19, 
2006, now Pat. No. 8,189,877. 

(60) Provisional application No. 60/728,733, filed on Oct. 
21,2005. 


(51) Int.Cl. 

G06K9/00 (2006.01) 

G06T 7/00 (2006.01) 

(52) U.S. Cl. 

CPC G06T 7/00 (2013.01); G06K 9/00657 

(2013.01) 

USPC 382/110 


(58) Field of Classification Search 

None 

See application file for complete search history. 

(56) References Cited 

U.S. PATENT DOCUMENTS 

5,565,998 A 10/1996 Coombs etal. 

6,484,099 B1 11/2002 Holzer-Popp et al. 

2002/0012071 Al* 1/2002 Sun 348/578 

2005/0104771 Al 5/2005 Terry etal. 

OTHER PUBLICATIONS 

Jiaguo Qi ; Cuizhen Wang ; Matricardi, E. ; Skole, D., “Improved 
selective logging detection with Landsat images in tropical regions”. 
Nov. 7, 2002, Geoscience and Remote Sensing Symposium, 2002. 
IGARSS ’02. 2002 IEEE International, vol. 4, pp. 2078-2080. 

* cited by examiner 

Primary Examiner — Atiba O Fitzpatrick 

(74) Attorney Agent, or Firm — Morgan Lewis & Bockius 

LLP 

(57) ABSTRACT 

The present invention provides systems and methods to auto- 
matically analyze Landsat satellite data of forests. The 
present invention can easily be used to monitor any type of 
forest disturbance such as from selective logging, agriculture, 
cattle ranching, natural hazards (fire, wind events, storms), 
etc. The present invention provides a large-scale, high-reso- 
lution. automated remote sensing analysis of such distur- 
bances. 

17 Claims, 10 Drawing Sheets 



0 


500 


1,000 


mv Km 

2,000 








Z Old 



CRITERIA 





















MONTE CARLO ANALYSIS 


U.S. Patent 


Mar. 31, 2015 


Sheet 4 of 10 


US 8,995,724 B2 






U.S. Patent 


Mar. 31, 2 


LU 

O 


0.40 


0.30 


6 0.20 


LU 

I 

Ll_ 

LU 


0.10 
0.00 


(A 


J 


0.5 


0.60 (B 
0.50 


o 

z: 

& 

o 


cc 


0.40 

0.30 

0.20 

0.10 

0.00 


0.5 


0.50 

0.40 


K~ 


o 


0.30 f 


y 0.20 


^ 0.10 
0.00 


0.5 










AUTOMCU 03-AUG-2001 


U.S. Patent 


Mar. 31, 2015 Sheet 7 of 10 


US 8,995,724 B2 






U.S. Patent 


Mar. 31, 2015 


Sheet 8 of 10 


US 8,995,724 B2 



FIG. 8 





U.S. Patent 


Mar. 31, 2015 


Sheet 9 of 10 


US 8,995,724 B2 



o CD 


o 

I — 
ZD 
< 
oes 

CD 

z 

CD 

CD 

O 

i 

CO 

< 

o 


DC CO 

dc 


CD LU 7? 


CD > a 

° 8 g 


LU LU 
O DC 
LU O 
DC LL 


O 

O 


o 

CO 


LU 

CD 



E 


CO 

LU 

CO 

CO 

< 

o 

CO 

LU 

Q 

O 

DC 

CL 


CO 

LU 

DC 

O 


o 

£ 

CO 

LU 

DC 

O 

LL 

LU 

□ 


CO 


CO 

c\j 







US 8,995,724 B2 


1 2 

REMOTE SENSING ANALYSIS OF FOREST BRIEF DESCRIPTION OF THE DRAWINGS 

DISTURBANCES 

For the purpose of illustrating the invention, there are 
STATEMENT REGARDING FEDERALLY depicted in the drawings certain embodiments of the inven- 

SPONSORED RESEARCH OR DEVELOPMENT 5 tion. However, the invention is not limited to the precise 

arrangements and instrumentalities of the embodiments 
This invention was supported in part by funds obtained depicted in the drawings, 
from NASA’s Large-Scale Biosphere-Atmosphere Experi- FIG. 1 depicts spatial distribution of selective logging in 
ment in Amazonia (LBA-ECO), grant number NCC5-675 five timber production states of the Brazilian Amazon for the 
(LC-21). The U.S. Government may therefore have certain 10 year intervals 1999-2000 (red), 2000-2001 (blue), and 2001- 
rights in the invention. 2002 (green). The states of Amazonas (AM), Amapa (AP), 

Tocantins (TO), Maranhao (MA), and the southern non-for- 
BACKGROUND OF THE INVENTION ested part of Mato Grosso were not included in the analysis. 

j - Light gray areas show the extent of indigenous reserves; dark 
Tropical forests have been threatened by increasing rates of ' gra y areas delineate federal conservation lands as of 1999 
deforestation or clear-cutting during the past three or more Instituto-Socioambiental. (Sao Paulo, Brazil, 1999) Map of 
decades (E. F. Lambin, H. J. Geist, E. Lepers, Ann. Rev. forest types, land-use change and protected areas in the Ama- 

Environ. Res. 28, 205 (2003)). Although deforestation, zon). 

largely for conversion of land to food crops or pastures, is the 20 FIG. 2 depicts a high resolution example of selective log- 
major destructive force in tropical forests worldwide, other g j ng results in 2001 -2002 from the CLAS processing in corn- 

forest disturbances such as the selective harvest of timber parison to deforestation mapping provided by the Brazilian 

have increased in frequency and extent (D. C. Nepstad et al., National Institute for Space Research (INPE (Instituto Nacio- 

Nature 398, 505 (1999), L. M. Curran et al., Science 303, na ] de Pesquisas Espaciais), “PRODES: Assessment of 

1000(2004)). In selective logging, a limited number of mar- 25 Deforestation in Brazilian Amazonia (http://www.obt.in- 
ketable tree species are cut, and logs are transported off-site to pe.br/prodes/index.html)” (2005)). 

sawmills. Unlike deforestation that is readily observed from FIG. 3 depicts the Carnegie Landsat Analysis System 
satellites, selective logging in the Brazilian Amazon causes a (CLAS) processing stream. 

spatially diffuse thinning of large trees that is hard to monitor FIG. 4 depicts the AutoMCU sub-model within CLAS, 

using satellite observations. Selective logging causes wide- 30 showing that each satellite image pixel is a calibrated reflec- 
spread collateral damage to remaining trees, sub-canopy veg- tance spectrum that is deconvolved into constituent fractional 

etation and soils, with impacts on hydrological processes, covers of photosynthetic vegetation (PV), non-photosyn- 
erosion, fire, carbon storage, and plant and animal species. thetic vegetation (NPV), and soil. Spectral endmemberlibrar- 

There is surprisingly little known about the extent or ies developed from extensive field and hyperspectral satellite 

impacts of selective logging throughout the tropical forests of 35 studies (TropiSpec) (Asner et al., 2005) are in a probabilistic 
the world, including the Amazon Basin. A survey of sawmills Monte Carlo unmixing approach to derive the percentage 
in the Brazilian Amazon suggested that 9,000-15,000 kin 2 of cover of PV, NPV and soil within each image pixel, 

forest had been logged in 1996-97 (D. C. Nepstad et al., FIG. 5 depicts spectral endmember bundles used in the 
Nature 398, 505 (1999)). The large uncertainty in this AutoMCU step of CLAS (from FIG. 3), which are (A) Pho- 

reported area resulted from necessary assumptions of the 40 tosynthetic vegetation, (B) Non-photosynthetic vegetation, 
wood volume harvested per area of forest. Sawmill surveys and (C) soil. Adapted from Asner et al. (2004a). 

can, at best, provide only a general idea of where and how FIG. 6 depicts an example of deforestation and water body 

much logging occurs because most operators buy timber at masking using Landsat thermal band 6 and the AutoMCU 

the mill gate rather than harvesting the wood themselves. result for photosynthetic vegetation (PV). 

Objective, spatially-explicit reporting on selective logging 45 FIG. 7 depicts an example of logging detection using 
requires either labor-intensive field surveys in frontier and CLAS. AutoMCU results from one year are differenced 

often violently contested areas, or by remote detection and against those of the next year. A directional pattern recogni- 

monitoring approaches. Previous studies of small areas show tion algorithm then uses the PV-change image to locate prob- 

the need for high-resolution observations via satellite. More- able logging decks, skids, and roads. 

over, most of the traditional analysis techniques employed for 50 FIG. 8 depicts a geographic coverage of study, showing the 
localized selective logging studies have proven insufficient Brazilian Legal Amazon with Landsat 7 satellite footprints, 

for large-scale selective logging assessments. A detailed FIG. 9 depicts an example showing how the CLAS logging 

comparison of Landsat satellite observations against field product is unique from the PRODES deforestation products 

measurements of canopy damage following selective logging provided by the Brazilian Space Research Institute, 
proved that traditional analytical methods missed about 50% 55 FIG. 10 depicts a block diagram of the CLAS system, 
of the canopy damage caused by timber harvest operations 

(G. P. Asner, M. Keller, R. Pereira, J. Zweede, Rem. Sens. DETAILED DESCRIPTION OF THE INVENTION 

Environ. 80,483 (2002)). 

All publications and patent applications herein are incor- 
BRIEF SUMMARY OF THE INVENTION 60 porated by reference to the same extent as if each individual 

publication or patent application was specifically and indi- 
The present invention provides systems and methods for vidually indicated to be incorporated by reference, 
automatically analyzing Landsat satellite data of forests. The The following description includes information that may 

present invention can easily be used to monitor any type of be useful in understanding the present invention. It is not an 
forest disturbance, such as, but not limited to, logging, agri- 65 admission that any of the information provided herein is prior 
culture, cattle ranching, natural hazards (fire, wind events, art or relevant to the presently claimed inventions, or that any 

storms), etc. publication specifically or implicitly referenced is prior art. 
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Unless defined otherwise, all technical and scientific terms 
used herein have the same meaning as commonly understood 
by one of ordinary skill in the art to which this invention 
belongs. Although any methods and materials similar or 
equivalent to those described herein can be used in the prac- 
tice or testing of the present invention, the preferred methods 
and materials are described. 

The computational analysis of Landsat Enhanced The- 
matic Mapper Plus (ETM+) satellite data was advanced using 
the new Carnegie Landsat Analysis System (CLAS) to detect 
and quantify the amount of selective logging in the major 
timber production States of the Brazilian Amazon. The 
approach provides automated image analysis using atmo- 
spheric modeling, detection of forest canopy openings, sur- 
face debris, and bare soil exposed by forest disturbances, and 
pattern recognition techniques. As discussed in greater detail 
below, CLAS provides detailed measurements of forest 
canopy damage at a spatial resolution of 30x30 meters, and it 
does so over millions of square kilometers of forest. 

CLAS was applied to five states — Para, Mato Grosso, Ron- 
donia, Roraima, and Acre — that account for -90% of all 
deforestation in the Brazilian Amazon. The analysis was con- 
ducted on a time-series of Landsat ETM+ imagery from 1 999 
to 2002. Across the five timber producing Brazilian states, the 
annual extent of selective logging ranged from 12,135 to 
20,651 km 2 (FIG. 1). These logging results represent new 
forest damage not accounted for in deforestation studies. 
Each year, the overlap between the results and the Brazilian 
National Institute for Space Research (INPE) annual defor- 
estation maps was only 6% (±5%). Moreover, only 19% 
(±1 1%) of the total area logged in any given year was subse- 
quently deforested three years later. Selective logging thus 
adds 60-128% more forest area damage than has been 
reported for deforestation alone in the same study period 
(Table 1). Selective logging was concentrated in the states of 
Mato Grosso and Para, where logging areas exceeded or 
nearly matched deforestation areas. In other smaller states, 
selective logging increased forest damage area by 10-35% 
over reported deforestation rates (Table 1). 


TABLE 1 


Selective-logging rates from 1999-2002 in five major timber-producing 
states of the Brazilian Amazon, with comparison to the deforestation 
rates reported bv INPE ('2005'). 

State 

1999-2000 rates 
(km 2 vr -1 ) 

Defor- 

Logged ested 

2000-01 rates 
(km 2 yr‘) 

Logged 

2001-02 rates 
(km 2 yr -1 ) 

Logged 

Logged 

Deforested 

Logged 

Deforested 

Acre 

64 

547 

53 

419 

111 

727 

Mato 

13,843 

6,176 

7,912 

7,504 

7,267 

6,880 

Grosso J 







Para 

5,939 

6,671 

5,343 

5,237 

3,791 

8,697 

Rondonia 

773 

2,465 

923 

2,673 

946 

3,605 

Roraima 

32 

253 

55 

345 

20 

54 

Total 

20,651 

16,112 

14,286 

16,178 

12,135 

19,963 


JOnly the northern 58% of Mato Grosso containing forested lands was included in the 
analysis. 


Conservation units such as indigenous reserves, parks and 
national forests generally afforded protection against log- 
ging. Elowever, exceptions included areas in northern Mato 
Grosso, where up to 880, 291, and 50 km 2 of logging were 
measured each year in the Xingu, Aripuana, and Serra 
Morena indigenous reserves, respectively (FIG. 1). In the 
southern portion of Para state, major logging disturbances 
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were observed in the Menkragnoti and Kayapo indigenous 
reserves, with up to 261 and 198 km 2 detected each year 
between 1999 and 2002. Federal forest reserves of Acre, 
Gorotire (Para), and Juruena (Mato Grosso) were harvested 
5 for timber at rates of up to 23, 90, and 380 km 2 each year, 
respectively. 

Extensive field validation studies showed that the canopy 
damage detection within CLAS is precise and accurate, as set 
forth below in the Materials and Methods section. Field vali- 
to dation studies showed false-positive and false-negative detec- 
tion rates of only 5%. Uncertainty caused by errors in atmo- 
spheric correction of satellite data, cloud cover, 
annualization, automated logging area delineation and 
manual auditing were 0.7-12.8% individually. After combin- 
15 ing all known sources of error, the analysis suggests an overall 
absolute uncertainty of up to 14% in total logging area. 

Selective logging contributes substantially to gross carbon 
fluxes from the Brazilian Amazon. Forest damage results 
from CLAS were combined with field-based forest canopy 
20 gap fraction and roundwood extraction data to calculate the 
total wood extraction rates. In 2000, 2001 and 2002, round- 
wood production averaged 49.8, 29.8 and 26.6 million cubic 
meters, respectively. The mean annual harvest intensities 
were 26.6, 21 .7 and 21 .4 m 3 ha -1 , which were generally lower 
25 than those reported by sawmill owners in 1 996. Nepstad et al. 
(1999) interviewed sawmill operators to estimate harvest 
intensities of 1 9, 28 and 40 m 3 per hectare in 1996. The total 
volume harvested equates to 10-15 million metric tons of 
carbon removed. In addition to roundwood, residual stumps, 
30 branches, foliage and roots are left to decompose in the forest, 
subsequently returning to the atmosphere as carbon dioxide 
over about a decade. The calculated average harvest intensity 
of 23.2 m 3 ha -1 equates to -8 Mg C ha -1 contained in round- 
wood, with an associated 34-50 Mg C ha' 1 of fine and coarse 
35 debris. The conversion of roundwood to carbon assumes an 
average wood specific gravity of 0.7 Mg m“ 3 and a propor- 
tional carbon content of 0.5 as in Keller et al. (2004a). Fallen 
debris creation was estimated based on data from Keller et al. 
(2004b) based on mean debris amounts found in logged for- 
40 ests (-30 nr ha' 1 harvested) subtracting the woody debris 
found in undisturbed forests. Upper and lower estimates were 
based on mean debris amounts plus root mean squared (RMS) 
error accounting for the uncertainty of estimates for both 
background and logged sites. Total debris was estimated as 
45 1 .4 times fallen debris to account for standing dead and roots 
(Keller et al., 2001). Integrated to the regional scale, the 
processing of roundwood and decomposition of residues lead 
ultimately to a gross flux of carbon from the forest of up to 
0.08 billion metric tons for each year of logging. The regional 
50 gross flux of carbon was estimated by multiplication of the 
range of carbon densities of debris created by the area logged. 
The range includes both variation in the annual area logged 
and uncertainty in the amount of debris created during log- 
ging. This value increases the estimated gross annual anthro- 
55 pogenic flux of carbon from Amazon forests by up to 25% 
over carbon losses from deforestation alone. Post-harvest 
forest regeneration reduces the net flux of carbon to the atmo- 
sphere below these values but the pace of regeneration after 
logging varies considerably. 

60 Selective logging doubles previous estimates of the total 
amount of forest degraded by human activities (Table 1), a 
result with potentially far-reaching implications for the ecol- 
ogy of the Amazon forest and the sustainability of the human 
enterprise in the region. In the future, improved monitoring of 
65 tropical forests will require high performance satellite obser- 
vations and new computational techniques. The results, pre- 
sented with explicit uncertainty analysis and transparency of 



US 8,995,724 B2 


5 

method, have located and quantified ubiquitous but previ- 
ously cryptic disturbances caused by selective logging. 
Definitions. 

As used herein, each of the following terms has the mean- 
ing associated with it in this section. 

The articles “a” and “an” are used herein to refer to one or 
to more than one (i.e. to at least one) of the grammatical object 
of the article. By way of example, “an element” means one 
element or more than one element. 

As used herein, the term “deforestation” refers to clear- 
cutting and conversion of the forest to other land uses, such as 
cattle pasture, crop agriculture, and urban and suburban areas. 

Without further description, it is believed that one of ordi- 
nary skill in the art can, using the preceding description and 
the following illustrative examples, make and utilize the com- 
pounds of the present invention and practice the claimed 
methods. The following working examples therefore, specifi- 
cally point out the preferred embodiments of the present 
invention, and are not to be construed as limiting in any way 
the remainder of the disclosure. 

EXPERIMENTAL EXAMPLES 

The invention is now described with reference to the fol- 
lowing examples. These examples are provided for the pur- 
pose of illustration only and the invention should in no way be 
construed as being limited to these examples but rather should 
be construed to encompass any and all variations which 
become evident as a result of the teaching provided herein. 

The materials and methods used in the experiments pre- 
sented in this Example are now described. 

Materials and Methods 
Processing Methodology 

The Carnegie Landsat Analysis System (CLAS) includes a 
general purpose computer programmed to use high spatial 
resolution satellite data for regional and global studies of 
forest disturbance. The computer system used is a multi- 
processor Linux system, but other systems can be used. 
CLAS is an automated processing system that includes: (i) 
atmospheric correction of satellite data; (ii) deconvolution of 
spectral signatures into sub-pixel fractional cover of live for- 
est canopy, forest debris and bare substrates; (iii) cloud, 
water, anddeforestationmasking; and (iv) pattern recognition 
algorithms for forest disturbance mapping. The following 
sections provide a detailed description of CLAS, illustrated 
by FIG. 3. 

Image Preparation and Atmospheric Correction 

The version of CLAS presented here ingests raw Landsat 
Enhanced Thematic Mapper Plus (ETM+) satellite imagery 
and applies sensor gains and offsets to convert from digital 
number (DN) to exo-atmospheric radiance. The radiance data 
are passed to a Hilly automated version of the 6S atmospheric 
radiative transfer model (Vennote et al.). The 6S program is 
integrated into the CLAS processing stream and uses monthly 
averages of aerosol optical thickness (AOT) and water vapor 
(WV) values from the Moderate Resolution Imaging Spec- 
trometer (MODIS) sensor onboard the NASA Terra space- 
craft. Time-stamping of MODIS AOT and WV data with 
Landsat data is done on an automated basis (FIG. 3). 
Sub-Pixel Analysis 

The CLAS process relies upon the quantitative determina- 
tion of fractional material cover at the sub-pixel scale (e.g., 
within each Landsat 30x30 m pixel). This core step employs 
a probabilistic spectral mixture sub-model that is run using 
the formulation shown in FIG. 4. This process spectrally 
decomposes each image pixel into fractional cover estimates 
(0-100% cover) of photosynthetic vegetation (PV) canopy. 
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non-photosynthetic vegetation (NPV), and bare substrate. 
This sub-model is based on an algorithm developed for forest, 
savanna, woodland and shrubland ecosystems. It is fully auto- 
mated and uses a Monte Carlo Unmixing (AutoMCU) 
5 approach to derive uncertainty estimates of the sub-pixel 
cover fraction values. The method uses three spectral end- 
member “bundles”, derived from extensive field databases 
and satellite imagery, to decompose each image pixel using 
the following linear equation: 
to 

P(^W=2/C , e -p(X.)J+e=/C i ,,.-p(X) pv +C„ ?v -p(U„ ; , v + 

^ substrate PPP substratell^^ 0) 

where p(X) e is the reflectance of each land-cover endmember 
(e) at wavelength k and 6 is an error term. Solving for the 
15 sub-pixel cover fractions (C e ) requires that the observations 
(pQ'-)pixei — in this case, Landsat ETM+ reflectance) contain 
sufficient spectral information to solve a set of linear equa- 
tions, each of the form in equation (1 ) but at different wave- 
lengths (L). 

20 Until recently, there were a limited number of spectral 
signatures of green and senescent vegetation and bare sub- 
strates for tropical regions. The mixture modeling technique 
requires spectral reflectance bundles (p^ v (X), p„ p Jk), and 
P substrateQ^)) that encompass the common variation in canopy 
25 and soil properties. Asner (1998) and Asner et al. (2003a, 
2004a) collected these spectral data using full optical range 
field spectroradiometers (Analytical Spectral Devices, Inc., 
Boulder, Colo., USA) during field campaigns conducted from 
1996 to 2000. The spectral endmember database encom- 
30 passes the common variation in materials found throughout 
the Brazilian Amazon, with statistical variability well defined 
(2004a). The bare substrate spectra have been collected 
across a diverse range of soil types, surface organic matter 
levels, and moisture conditions. Spectral collections for NPV 
35 have included surface litter, senescent grasslands, and defor- 
estation residues (slash) from a wide range of species and 
decomposition stages. 

In contrast to the NPV and bare substrate spectra that can 
be collected via ground-based spectroscopic measurements, 
40 the photosynthetic vegetation (PV) spectra of forest species 
require overhead viewing conditions. This is very difficult in 
forest canopies with heights typically ranging from 10-50 m. 
Spectral measurements of individual leaves, stacks of foliage, 
or partial canopies (e.g., branches) introduce major errors in 
45 spectral mixture models and cannot be used (Asner, 1998). 
Therefore, canopy spectra were collected using the Earth 
Observing- 1 (EO-1) Hyperion sensor, the first spaceborne 
hyperspectral sensor for environmental applications (Ungar 
et al.). The PV spectral bundle was derived from more than 
50 40,000 spectral observations made at 30 m spatial resolution 
with Hyperion (images taken throughout 1999), atmospheri- 
cally corrected to apparent top-of-canopy reflectance using 
the ACORN-4 atmospheric correction algorithm for hyper- 
spectral data (ImSpec Inc., Palmdale, Calif. USA), and con- 
55 volved to Landsat ETM+ optical channels (Asner et al., 
2005). These green vegetation spectra thus inherently 
included the variable effects of intra- and inter-crown shad- 
owing, which are prevalent in tropical forests (Gastellu- 
Etchegorry et al.). In Amazonia, shade fractions average 25% 
60 cover in humid tropical forests, but the variance is high with 
standard deviations of 12% or more (Asner et al., 2003b). 

It is thus critically important to note that the PV results 
include shade, which varies substantially with forest struc- 
ture. Using a separate shade endmember is attractive (Souza 
65 et al., 2000), but doing so with multi-spectral Landsat data 
and such high shadow fraction variability often results in an 
under-determined spectral and mathematical problem in lin- 
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ear mixture models . That is, there are many viable solutions to 
the mixture modeling problem in forests. Imaging spectros- 
copy (hyperspectral) data are needed to solve this problem 
(Roberts et ah, 1 993). This issue was avoided by accepting the 
limitations of incorporating variable shade directly into the 5 
PV bundle derived from the EO-1 Hyperion sampling of 
undisturbed forest canopies in Brazil. The PV bundle 
includes spectra from mature forest, late-stage forest 
regrowth, and logged forest of at least five years post-harvest. 

In the end, the total number of spectra retained in the end- 10 
member bundles fortheAutoMCU sub-model was 252, 611, 
and 434 for PV, NPV and bare substrate, respectively (FIG. 

5). These spectra represent more than 130,000 field and spa- 
cebome spectrometer observations collected over a five-year ( 5 
period of study (Asner et ah, 2005). 

Non-Forest Masking and Atmospheric Compensation 
A series of automated masks were designed to exclude 
clouds, water bodies, cloud shadows, non-image and non- 
forest areas (e.g., pasture, urban and agriculture) from the 20 
CLAS processing stream (FIG. 3 ). Prior to execution of the 
AutoMCU sub-model, clouds are masked using the thermal 
channel (band 6) from the raw Landsat images. Asner et al. 
(2005) found that a thermal band threshold DN value of 125 
can conservatively detect cloudy pixels over Amazonia. 25 
Water bodies are masked by finding pixels in the calibrated 
Landsat reflectance data in which bands 1 -4 (blue, green, red, 
and near-infrared) have a negative slope. Only water displays 
such a negative reflectance slope with increasing wavelength. 
Non-image areas containing zero values are also masked. 30 
Cloud shadows are identified using the root mean square 
error (RMSE) image that results from the AutoMCU process- 
ing (FIGS. 3 - 4 ). Areas shadowed by clouds have large RMSE 
values and are masked by identifying pixels above a specific 
RMSE threshold (Asner et al., 2005). To limit the logging 35 
analysis to forested areas, Landsat thermal band 6, combined 
with the AutoMCU results, is used to identify pixels contain- 
ing primarily forest and non-forest areas. Forests have a lower 
brightness temperature and a higher PV fractional cover than 
deforested lands. A conservative PV fractional cover thresh- 40 
old of 60% was employed to delineate forest cover in the PV 
mask. The minimum and maximum thermal thresholds, 
which encompass forested areas in the thermal mask, are 
dynamically generated for each image by calculating the 
mean thermal value of all pixels having a PV fraction cover 45 
greater than 80% and then masking all pixels with values >15 
digital numbers (DN) from the mean thennal value. These 
final masking steps have the added feature of removing 
residual clouds and cloud shadows that were missed in the 
masks applied earlier in the CLAS process (FIG. 6). 50 

Although atmospheric correction was performed on the 
raw imagery before processing through the AutoMCU sub- 
model, residual atmospheric effects can persist (Asner et al., 
2005). These residual effects exist spatially within a scene 
and temporally between scenes. These effects were greatly 55 
reduced prior to automated logging detection (next section) 
by calculating the average change in fractional forest cover in 
55 km 2 subsets of the imagery. These large geographic sub- 
sets are made at a spatial scale far greater than that of the most 
extensive logging activities, so temporal differences in the 60 
overall forest fractional cover at this scale are a result of 
atmospheric effects (e.g., haze) or forest phenology. These 
false fractional cover changes are normalized by adjusting the 
background forest temporal variation to zero. Since distur- 
bances related to logging or other anthropogenic activities 65 
occur at a much smaller spatial scale than is considered in this 
processing step, normalization of the forest values across 


8 

large areas does not affect the CLAS process in discriminat- 
ing true disturbances from the surrounding forested areas. 
Pattern Recognition 

The specific criteria used in this procedure were deter- 
mined following a comprehensive analysis and review of the 
forest responses to logging at various intensities in the Bra- 
zilian states of Para, Mato Grosso and Acre where field stud- 
ies were conducted. The mean and standard deviation frac- 
tional cover images from the AutoMCU step in CLAS 
provide quantitative data on canopy damage and forest dis- 
Uirbance intensity from which selectively logged areas can be 
determined (FIG. 3). By identifying areas of canopy distur- 
bance that are arranged in specific spatial patterns, it is pos- 
sible to detect logged areas on an automated basis. The pri- 
mary method by which logging is detected is image 
differencing, where pairs of AutoMCU sub-pixel fractional 
cover images, separated by approximately one year, are used 
to create images of PV (forest canopy) and NPV (surface 
woody and senescent vegetation material) change that indi- 
cate areas of relative canopy disturbance or recovery. Forest 
disturbances in these images always have reductions in PV, 
simultaneous with increases in NPV fractional cover. 

Logging activity results in low intensity forest disturbances 
from tree felling gaps, moderate intensity linear features from 
skid trails along which felled trees are dragged by tractors or 
skidders, and high intensity points of damage called log decks 
where logs are loaded onto trucks for transportation. The log 
decks are connected by logging roads, seen as linear features 
causing large reductions in the fractional cover of PV, to local 
roads or rivers for transportation to markets. These patterns 
are unique to logging throughout most of the Amazon, and 
thus they serve as the basis upon which the method for log- 
ging detection functions. CLAS identifies points (e.g., treefall 
gaps and log decks) and linear features (e.g., skid trails and 
logging roads) of recent disturbance occurring in forested 
areas. As these features also exist at a lower frequency in 
intact forest regions, their spatial density and diversity (see 
definition in next section) are calculated to identify those 
areas having disturbances in patterns most indicative of log- 
ging activity. The procedure then identifies these areas for 
further analysis by creating point maps, termed logging 
nodes, indicating their locations. 

Log decks are automatically detected by searching for 
pixels where PV decreases significantly in a 30 m pixel cen- 
tered on a 7x7 pixel kernel (4.41 ha). A positive detection is 
flagged when pixels with large PV reduction are surrounded 
by three concentric rings of incrementally greater PV cover 
surrounding the target pixel. This indicates an increase in 
canopy damage with greater proximity to the log deck, a 
pattern consistent with most logging activities. 

The strategy for detecting decks works well in areas logged 
at higher intensities, as the decks tend to be abundant and 
equally spaced. However, in areas where the logging is more 
haphazard, where the forest damage is extremely high or low, 
or where the roads themselves also function as loading zones, 
individual log decks are not always distinguishable. Skids 
trails are a typological feature of selective logging practices, 
and they are the single-most ubiquitous surface feature found 
in harvested areas (Pereira et al., 2002; Asner et al., 2004). 
The presence of skid trails is quantifiable based on large 
decreases in PV fractional cover in linear or near-linear pat- 
terns (Asner et al., 2004a). To detect the concentration of skid 
trails and auxiliary roads, a moving 6x6 pixel (3 .24 ha) kernel 
is applied to the PV change image to enhance linear features 
in the N-S, E-W, NE-SW, and NW-SE directions (FIG. 3). 
The number of directions in which the linear features are 
arranged (which are defined herein as their diversity), and 
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their spatial density, in conjunction with the presence or 
absence of logging decks, is calculated for each location. 
With this information, it is possible to automatically distin- 
guish probable logging events. In general, areas of greater 
logging intensity have a roughly equal proportion and higher 5 
density of linear features with the presence of logging decks. 
Lower intensity areas are normally dominated by one direc- 
tion of linear feature and have few or no logging decks. An 
example of a typical logging detection is shown in FIG. 7. 

Final Integration to 

After the linear and logging deck pattern recognition steps 
are completed, CLAS automatically integrates the various 
results to identify contiguous pixel clusters of probable log- 
ging activity. This process starts by creating a list of the 
logging nodes that are identified in the previous steps . Logged 1 5 
areas are identified using a moving kernel approach. A base 
kernel of 7x7 pixels (4.41 ha) and four 3x3 pixel (0.81 ha) 
subset kernels, one located at each comer of the base kernel, 
are used. The base kernel begins at each logging node and 
tests the criteria described below. If the area in question tests 20 
positive, the analysis kernel is moved to its 7x7 pixel neigh- 
bors to the north, south, east, and west, which are then each 
tested against the criteria (FIG. 3). This iterative process 
continues until all neighbors have been evaluated or the maxi- 
mum logged cluster size (maximum of 17 positive detections 25 
per logging node) has been reached. The input layers and 
specific criteria tested within the base and subset kernels are 
described below. For the criteria below, all units for PV and 
NP V are % fractional cover within a pixel ; units for PV Cl and 
NP V Cl are % change in cover fractions between image dates . 30 
Input Layers to Logged Area Detection Procedure: 

Logging node map 

Thermal RMS mask (dynamically generated in earlier pro- 
cedure) (T-mask) 

PV mask (>60% fractional cover) (PV-mask) 35 

PV change difference image (PV Cl) 

NPV change difference image (NPV Cl) 

After image PV (AI PV) 

Base Kernel Criteria: 

75% good data pixels (not cloud, cloud shadow, or water) 40 
Non-forested area <0.54 ha (12.2%); based on T- and PV- 
masks. 

60%<Mean AI PV>93% 

Mean PV CI>-9% 

Mean NPV CI<2% 45 

Mean PV Cl standard deviation>33% 

Mean NPV Cl standard deviation>46% 

More than 6 pixels (0.54 ha) with PV Cl values>80% 

More than 6 pixels (0.54 ha) with NPV Cl values<-85% 
Masked area<0. 1 8 ha 50 

Subset Kernel Criteria: 

>2 subsets with PV CI>32% standard deviation 
>2 subsets with mean PV CI<3% and >60% 

>2 subsets having >1 pixel (0.09 ha) with a PV Cl 
value>80% 55 

>2 subsets with NPV CI>46% standard deviation 
>2 subsets with mean NPV CI<-5% and >-65% 

>2 subsets having >1 pixel (0.09 ha) with a NPV Cl 
value<-85% 

Manual Audit 60 

Maps of probable logging events were visually audited to 
verify whether an area is being logged or not, in accordance 
with criteria established for identification of logged areas (see 
criteria below). In this process, false positives and negatives 
were manually removed and added. In this Amazon study, 65 
two analysts were employed during the audit, and their results 
and uncertainties were monitored and compared. 
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The audit logging criteria are divided into high- and low- 
damage obvious and non-obvious categories. These catego- 
ries encompassed all probable logging events in the study 
area and were identified after extensive review of logging 
events identified in the field. The criteria applied in each 
category are listed below. 

High-Damage Obvious Criteria 
Abundance of logging decks 

Obvious linear features (including primary -tertiary access 
roads and skid trails) 

Severe canopy damage visible in PV change difference 
image 

Areal extent normally >1 ha 

Evidence of logging from previous years in close proxim- 
ity (<-15% PV change difference image) 
High-Damage Non-Obvious Criteria 
Few to no logging decks 
Few to no linear features 

Severe canopy damage visible in PV change difference 
image (>15%) 

Presence of access roads or rivers, if not adjacent to a 
anthropogenic non- forest areas 
Areal extent normally >1 ha 

Evidence of logging from previous years in close proxim- 
ity (<-15% in PV change difference image) 
Low-Damage Obvious Criteria 
Few to no logging decks 
Obvious linear features 
Presence of access roads or rivers 
Often tree-like in formation (graduating from higher to 
lower damage linear features) 

Not a linear feature connecting non-forest areas, or other- 
wise used for general transportation 
Often encompassing large extents 
Evidence of logging from previous years in close proxim- 
ity (<-15% PV change difference image) 

Low-Damage Non-Obvious Criteria 
Few to no logging decks 
Few to no linear features 

Close proximity to access (i.e. roads, rivers or anthropo- 
genic non- forest areas) 

Speckles of recent canopy damage (felling gaps; >15%) in 
PV change difference image occurring at a density 
greater than in the surrounding forest areas 
Areal extent normally <6.5 ha 

Evidence of logging from previous years in close proxim- 
ity (<-15% in PV change difference image) 

Variation in the final logging products may result from 
differences in observer application of the manual audit pro- 
cedure. The error associated with these user-specific differ- 
ences was tested by randomly selecting 25 subsets (400x400 
pixels=144 km 2 ) of the logging results throughout the Ama- 
zon study area. These areas were then manually audited by 
two separate observers — one with in-depth experience with 
the auditing procedure, and one with no experience whatso- 
ever. Paired t-tests were then used to assess differences in the 
final auditing results. This test showed no statistical differ- 
ence between auditors (p=0.21). As a result, it is concluded 
that the auditing step used herein — the only significant non- 
automated step in the CLAS processing stream — is not sen- 
sitive to user biases. 

Annualization 

To determine the total amount of logging within a Brazilian 
state for a given year interval, it is necessary to account for the 
seasonality of logging operations during the year. Since the 
Landsat ETM+images are selected based on lowest possible 
cloud cover, the time interval for which the logging activity is 
identified in the pair of images is not exactly one year, and so 
may contain more or less than one year of logging. Therefore, 
a time period was defined to which annual logging estimates 
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must be standardized. The annualization methodology used 
by the Monitoring the Brazilian Amazon Gross Deforestation 
Project (PRODES) was adopted. Assuming that logging, like 
deforestation, occurs during the dry season, the annual log- 
ging rate was determined by prorating the amount of logged 
area to a complete dry season. This standardization depends 
upon the onset and length of the dry season for each image, 
based the geographic location of the scene center. The general 
spatial pattern of the Amazon dry season was provided by 
Marengo et al. (2001 ), based on the daily amount of outgoing 
longwave radiation. 

Geographic Coverage 

The Brazilian Amazon basin covers an area of approxi- 
mately 4.1 million km 2 . Analysis of the entire region with 
Landsat Enhanced Thematic Mapper-Plus (ETM+) imagery 
would require approximately 220 scenes per year or 880 
images for the years 1 999-2002, yet much of the northwestern 
Amazon still contains relatively little deforestation and log- 
ging (Nepstad et al., 1999). Therefore, the study was limited 
to the States of Acre, Para, Mato Grosso (northern 58% of the 
state containing most of the forested area), Rondonia, and 
Roraima (FIG. 8). These five states contain -90% of the 
deforestation reported by Brazil for all of the Legal Amazon, 
and thus are the most important areas for logging studies 
today. This strategy reduced the number of required Landsat 
images to 480 scenes. 

The geographic overlap was evaluated between the Brazil- 
ian Space Research Institute’s deforestation maps and the 
CLAS logging maps for all Landsat images used in the study. 
An example of this comparison is shown in FIG. 9. It was 
found that the logging detections overlapped with PRODES 
deforestation maps only 6% (±5%) of the time in any given 
year. Up to three years following harvest, a maximum of 19% 
(±1 1%) of logged areas were subsequently deforested (clear 
cut). Therefore, these results are not redundant with defores- 
tation and thus represent forest damage that has been unac- 
counted for in previous State- and Basin-scale forest distur- 
bance estimates. 

The study covered the period 1999 to 2002, which is prior 
to the failure of the Scan Line Corrector (SLC) in the ETM+ 
instrument onboard Landsat 7. Following the SLC failure, 
roughly 40% of each acquired Landsat image is missing data. 
To seek out alternatives to Landsat 7, a satellite inter-com- 
parison of logging detection capability based on the network 
of low- and high-intensity logging sites in Amazonia was 
conducted (Asner et al., 2004). Comparisons were made 
among the detection capabilities of hyperspectral (EO-1 
Hyperion), multi-view angle (Terra-MISR), high spatial reso- 
lution multi-spectral (EO-1 Advanced Land Imager, Landsat 
5 Thematic Mapper, Landsat 7 ETM+, CBERS-2, SPOT, 
Terra- ASTER- VNIR), and low spatial resolution multi-spec- 
tral (Terra-MODIS, AVHRR) data. The only sensor to meet or 
exceed the performance of Landsat ETM+was EO-1 Hype- 
rion; all others failed to detect at least 80% of the logging 
damage in the field sites. However, EO-1 Hyperion is a hyper- 
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spectral technology demonstration with extremely limited 
spatial and temporal coverage, making its application to 
large-area analysis intractable. It was concluded that the com- 
bination of spectral resolution (6 optical channels), spectral 
5 signal-to-noise perfonnance, and spatial resolution was criti- 
cally important in detennining the amount of logging that 
could be reliably detected. Landsat 7 ETM+ provides the 
minimum performance needed for reliable analysis of selec- 
tive logging in Amazon forests, and improved space-based 
10 technology is critically needed to remedy the current limita- 
tions. 

Unobserved Areas 

In some cases, there is a lack of sufficiently cloud-free 
j - imagery to determine logging extent using the pattern recog- 
nition portion of CLAS. For example, cloud cover precluded 
the use of two image pairs (Landsat path-row 222/062 and 
226/062) in the state of Para, one in 1999 and the other in 
2000. To estimate logging in these areas, it was necessary to 
20 employ the single-scene analysis approach detailed by Asner 
et al. (2004a), which demonstrated that manual interpretation 
of single-date AutoMCU results have a temporal sensitivity 
of about one year following logging. Therefore, the extent of 
new logging in these two Paraimages as a proxy for annual 
25 logging rates in these regions was estimated. 

In other cases, persistent cloud cover resulted in no images 
for certain areas for certain years. In these instances, the 
amount of logging for the area was estimated from Landsat 
observations made in the closest year. Of the 480 Landsat 
30 images employed throughout the entire study, this was nec- 
essary 10.5% of the scenes. 

Uncertainty Analyses 

Any large-scale, complex remote sensing study must track 
and manage sources of uncertainty in the final results. This is 
35 important because there are many steps that can lead to errors. 
The uncertainty was carefully quantified in four key areas: (i) 
atmospheric correction (aerosol and water vapor), (ii) unob- 
served areas caused by persistent cloud cover, (iii) annualiza- 
tion. and (iv) auditor uncertainty. 

40 Atmospheric Uncertainty 

In the CLAS processing stream, Landsat ETM+ images are 
atmospherically corrected using the 6S atmospheric correc- 
tion algorithm (Vermote, 1997), with monthly averages of 
aerosol and water vapor inputs from the MODIS satellite 
45 sensor. The sensitivity of the CLAS AutoMCU algorithm to 
atmospheric correction errors was comprehensively assessed 
by Asner et al. (2005), and was found to be minimally sensi- 
tive to uncertainties in aerosol and water vapor from MODIS. 
To further understand the effect that the atmospheric correc- 
50 tion has on the sensitivity of entire CLAS process, five Land- 
sat image pairs were atmospherically corrected using ran- 
domly-selected, monthly aerosol and water vapor values 
from MODIS. The difference in the amount of automatically 
detected logging between the different atmospherically-cor- 
rected images was only 0.7% (Table 2). 

TABLE 2 


Primary sources of uncertainty in CLAS analyses of selective logging 
extent in forests of Amazonia. 


Percentage of 

Source Total Logged Area Method 


Atmosphere ±0.7% Difference in automatically detected logged 

area between the atmospherically corrected 
image and an image with randomly selected 
atmospheric characteristics. 
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TABLE 2 -continued 


Primary sources of uncertainty in CLAS analyses of selective logging 
extent in forests of Amazonia. 


Percentage of 


Source 

Total Logged Area 

Method 

Unobserved Area 

+5% 

Percentage of cloud- and shadow-covered area 
compared to total logged area 

Annualization 

±2-9% 

Standard error of the difference between dry 
season length for matched pairs of consecutive 
years from 1979-1996. 

Auditor 

±12.8% 

Standard error of difference between auditor 
estimates, on a per km 2 of logging basis 

TOTAL ESTIMATED 
ERROR 

±11-14% 

Root mean square error 


Unobserved Area Uncertainly 

When cloud and cloud-shadow cover is greater than 50% in 
any 5,625 kni2 area (2,500x2,500 pixels), the area of 20 
observed logging is used to estimate the amount of logging in 
the unobserved, cloudy areas. The sensitivity to this type of 
error was assessed by simply quantifying the fractional cover 
of clouds and cloud shadows in comparison to observed log- 
ging extent. The calculated absolute uncertainty caused by 22 
this step was approximately +5% over the five states (Table 
2 ). 

Annualization Uncertainty 

Although the rate of logging is assumed constant through- , f) 
out the dry season, there is a level of uncertainty inherent in 
this assumption. Marengo et al. reported rainy season length 
for five regions of the Amazon (i.e., North Amazonia, Central 
Amazonia, Mouth of Amazon, Southeast Amazonia, and 
Southwest Amazonia) for the period 1979-1996. To deter- 35 
mine the uncertainty in the logging estimate related to 
assumption of dry season length, a series of matched pairs of 
dry season length for two consecutive years (e.g., 1 979- 1 980, 
1980-1981, 1995-1996) was compiled to calculate the stan- 40 
dard error of the difference in dry season length for each 
region. This standard error (in days) was divided by the aver- 
age length of the dry season for the respective region to 
express the uncertainty in percent of dry season. This percent- ^ 
age uncertainty was then applied to actual satellite image 
pairs or by averaging the uncertainty for states located 
between two regions. These uncertainties ranged from 2-9% 
as a result of interannual variation in dry season length (Table 
2 ). 50 

To further assess the sensitivity of the logging area esti- 
mates to the annualization and timing of the dry season, the 
estimates were also annualized without the constraint that 
logging activity only occurs only during the dry season. These 55 
results are reported in Table 3, with comparison to the pre- 
ferred results that appear in the main text (Table 2). It is clear 
that the differences between these two assumptions can be 
large in the smaller states (e.g.. Acre, Rondonia), where the 
estimate of logged areas is more sensitive to the acquisition 60 
dates of a smaller number of annual satellite image pairs. 
However, in the larger states, these uncertainties tend to bal- 
ance out. In the majority of cases, the amounts of logging 
estimated without the dry season constraint still falls within 65 
the minimum and maximum limits of estimated logged area 
caused by other sources of uncertainty (discussed below). 


TABLE 3 


Logging estimates for Brazilian states in the Amazon using the 
dry-season annualization protocol from INPE (Marengo et al., 2001) 
and a calendar-year annualization. 



1999-2000 rates 
(knTvrb 

2000-2001 rates 
(km 2 vr - L 

2001-2002 rates 
(km 2 vr"b 

State 

Dry- 

season 

Protocol 

Calendar 

Year 

Dry- 

season 

Protocol 

Calendar 

Year 

Dry- 

season 

Protocol 

Cal- 

endar 

Year 

Acre 

64 

91 

53 

48 

111 

117 

Mato 

13,843 

11,762 

7,912 

7,783 

7,267 

7,182 

Grosso 

Para 

5,939 

5,030 

5,343 

5,159 

3,791 

3,751 

Rondonia 

773 

694 

923 

902 

946 

638 

Roraima 

32 

32 

55 

55 

20 

20 

TOTAL 

20,651 

17,609 

14,286 

13,947 

12,135 

11,708 


Auditor Uncertainty 

Each auditor reviewed a set of the same 25 image subsets 
(400 by 400 pixels) in which most images include some form 
of logging. A test was performed in which a novice and an 
experienced image analyst manually delineated areas con- 
taining logged forest. This comparison was used to calculate 
one standard error of the difference in logging assessments 
between auditors for each image subset. The standard error 
between auditors was 0.69 km 2 of logging, which when 
scaled by the average amount of logging identified by the two 
analysts (5.4 km 2 ), resulted in an uncertainty of 12.8% (Table 
2 ). 

These different sources of uncertainty were compiled and 
used to estimate an overall uncertainty in the logging extent 
estimates of 11-14% for each Brazilian state in each year of 
analysis (Table 2). These uncertainties were then propagated 
to the Basin scale for annual estimates of selective logging for 
the years 2000, 2001 , and 2002 (Table 4). 

TABLE 4 


Minimum-maximum logging estimates for Brazilian states in the 



1999-2000 rates 
(km 2 yr -1 ) 
Logged 

2000-01 rates 
(km 2 yr -1 ) 
Logged 

2001-02 rates 
(km 2 yr _1 ) 
Logged 

State 

Mini- 

mum* 

Maxi- 

mum! 

Mini- 

mum* 

Maxi- 

mum! 

Mini- 

mum* 

Maxi- 

mum! 

Acre 

Mato Gross o| 

54 

11,801 

78 

16,521 

45 

6,744 

66 

9,481 

94 

6,195 

133 

8,453 
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TABLE 4-continued TABLE 5-continued 


Minimum-maximum logging estimates for Brazilian states in the 
Amazon based on uncertainties in CLAS logging methodology. 


State 

1999-2000 rates 
(km 2 yr -1 ) 
Logged 

2000-01 rates 
(km 2 yr -1 ) 
Logged 

2001-02 rates 
(knEyr -1 ) 
Logged 

Mini- 

mum* 

Maxi- 

mum^ - 

Mini- 

mum* 

Maxi- 

mum^" 

Mini- 

mum* 

Maxi- 
mum j" 

Para 

4,905 

7,419 

4,421 

6,536 

3,138 

4,844 

Rondonia 

657 

931 

785 

1,076 

804 

1,113 

Roraima 

27 

38 

46 

66 

17 

26 

Total 

17,444 

24,987 

12,041 

17,225 

10,248 

14,569 


^Composed of atmospheric, temporal interpolation, annualization, and auditor uncertainties 
(see text for definitions). 

f Includes all uncertainties plus cloud interpolated area. 

J Includes Northern Mato Grosso only. 

Validation 

A comprehensive validation study of the logging extent 
results derived from the CLAS processing stream was carried 
out. Previous validation studies were highly detailed (tree- 
by-tree level) damage assessments, but were limited to fewer 
logging sites and only to the Auto MCU portion of the process 
(Asner et al., 2004a). Pereira et al., Asner et al. (2002b, 
2004b), and Keller et al. (2004b) carried out extensive field 
studies from 1997-2002 in logging areas subjected to a wide 
range of harvest methods, intensities, and canopy damage 
levels. These studies included the development of high-reso- 
lution global positioning system (GPS) coverages of logging 
extent in conventional and reduced-impact logging sites in 
eastern Para, central Para, and northern Mato Grosso. These 
areas contained the most intensive and widespread logging in 
the entire study. 

There were a total of 45 harvest/image combinations avail- 
able for this validation study. The images areas were only 
considered where the harvest blocks were free of clouds and 
whose harvest month was known when that knowledge was 
essential. All logging events were contained in three Landsat 
images: Fazenda Cauaxi in eastern Para (Landsat path/row 
223/063); Tapajos National Forest in central Para(path/row 
227/062); and Juruena in northern Mato Grosso (path/row 
229/067). The timber harvest dates of areas contained within 
these images ranged from 1 997 to 2002 (Table 5). About half 
of the logging sites were harvested using conventional (high- 
damage) techniques, and the other half employed reduced- 
impact (low-damage) logging methods (Asner et al., 2004b). 
Logging areas ranged in size from 1 1 to 1,079 ha. This wide 
range of logging block sizes and canopy damage levels pro- 
vided a substantial geographic data set against which to test 
CLAS. 

TABLE 5 


Validation of CLAS logging detection method. 

Logging Harvest Image Harvest % Logging 

Block Type Date Date Detected 

Harvests more than 12 months prior to satellite 
imaging or after satellite imaging 


Cauaxi 

CL 

2000 

1996 

0 

Cauaxi 

RIL 

2000 

1996 

0 

Cauaxi 

RIL 

2000 

1998 

0 

Cauaxi 

CL 

2000 

1998 

13 

Cauaxi 

CL 

2000 

2000 

14 

Juruena 

CL 

2000 

1998 

0 

Juruena 

CL 

2000 

1999 

5 


Validation of CLAS logging detection method. 


Logging 

Harvest 

Image 

Harvest 

% Logging 

Block 

Type 

Date 

Date 

Detected 

Juruena 

CL 

2000 

2001 

16 

Juruena 

RIL 

2000 

2002 

9 

Tapajos 

RIL 

2000 

1997 

0 

Tapajos 

RIL 

2000 

2002 

0 

Tapajos 

RIL 

2000 

2000 

8 

Jumena 

RIL 

2001 

2002 

1 

Tapajos 

RIL 

2001 

2002 

0 

Cauaxi 

CL 

2001 

1996 

0 

Cauaxi 

RIL 

2001 

1996 

0 

Cauaxi 

RIL 

2001 

1998 

0 

Cauaxi 

CL 

2001 

1999 

0 

Cauaxi 

RIL 

2001 

1999 

0 

Cauaxi 

CL 

2001 

1998 

4 

Juruena 

CL 

2001 

1998 

0 

Juruena 

CL 

2001 

1999 

38 

Tapajos 

RIL 

2001 

1997 

0 

Tapajos 

RIL 

2001 

1999 

0 

Tapajos 

RIL 

2001 

2000 

0 

Cauaxi 

CL 

2002 

1996 

0 

Cauaxi 

RIL 

2002 

1996 

0 

Cauaxi 

RIL 

2002 

1998 

0 

Cauaxi 

CL 

2002 

1999 

0 

Cauaxi 

RIL 

2002 

1999 

6 

Cauaxi 

RIL 

2002 

2000 

8 

Cauaxi 

CL 

2002 

1998 

0 

Cauaxi 

CL 

2002 

2000 

34 

Jumena 

CL 

2002 

1998 

0 

Jumena 

CL 

2002 

1999 

0 

Jumena 

RIL 

2002 

2000 

0 

Tapajos 

RIL 

2002 

1999 

0 

Tapajos 

RIL 

2002 

2000 

0 

Tapajos 

RIL 

2002 

2000 

16 

Harvests less than 1 2 months 

prior to satellite 

imaging 

Cauaxi 

CL 

2000 

1999 

75 

Cauaxi 

RIL 

2000 

1999 

59 

Tapajos 

RIL 

2000 

1999 

91 

Tapajos 

RIL 

2000 

2000 

84 

Cauaxi 

CL 

2001 

2000 

44 

Tapajos 

RIL 

2002 

2002 

36 


Geographic information system (GIS) coverages of the 
logging areas listed in Table 5 were overlaid on the CLAS 
products, and statistical data on logging detection percentage 
area and logging type were calculated. Results were orga- 
nized by success or failure in detecting the logged areas and 
their spatial extent. For analysis purposes, a false-negative 
detection was declared when CLAS missed areas logged in 
the 1 2 months prior to image date. A false-positive detection 
was declared for areas not logged in the 12 months prior to 
image date and when more than 25% of the false area was 
detected as logging. The 12-month limit, a result of the acqui- 
sition dates of the before and after Auto MCU images used in 
the change differencing process, was selected based on the 
known sensitivity of the AutoMCU algorithm within CLAS 
(S3). 

Of the 45 image combinations tested, only two false-posi- 
tives and two false-negatives occurred (Table 5). For the two 
false-positives, the detected logging areas ranged from 
34-38% of the true area and were logging blocks that had been 
harvested two years prior to satellite imaging. Therefore, it is 
considered that these two blocks as false-positives only in the 
sense that it was not intended for CLAS to find logging sites 
that are more than a year old. Further review of the false- 
positive from Cauaxi pointed to issues of geo -registration of 
exact harvest boundaries between the GIS and the imagery. 
The other false-positive (Juruena) was very clearly re-har- 
vested over roughly 40-50% of the originally delineated log- 
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ging block (which had been harvested three years prior to 
re-harvest), hence it was not considered further. For the two 
false-negatives, the detected logging areas ranged from 
36-44% of the true area. Further review of the false-negative 
from Cauaxi also revealed co-registration problems between 
field data and remote sensing imagery. Hie other false-nega- 
tive (Tapajos) occurred in a large block where the harvesting 
took place for several months, which caused the first har- 
vested portions to be too old (14 months) for detection due to 
regrowth. 

A chi-square test of the results in Table 4 shows that there 
is a highly significant association between logging and CLAS 
algorithm detection (x 2 =17.0>10.8; p<0.001 (Table 6) This 
test is conservative because it is mainly testing the sensitivity 
to currently versus previously logged areas, and not to intact 
forest that have not been harvested. 

TABLE 6 

Chi square (yf) test of the validation results 
of the CLAS logging methodology. 

Logged Logged 

No Yes 

Observed 


Detected 

No 

37 

2 

Detected 

Yes 

Expected 

2 

4 

Detected 

No 

33.8 

5.2 

Detected 

Yes 

5.2 

0.8 


* Number of harvest/image combinations that met the following criteria: harvest months 
were known when that knowledge was essential for validation; detection was not obviously 
interfered with because of proximity of clouds/cloud shadows. 

Roundwood Estimates from Remote Sensing and Field Stud- 
ies 

The roundwood production was calculated from logging 
areas detected using the CLAS approach, combined with 
field-based relationships between remotely sensed canopy 
opening (PV cover), forest canopy gap fraction and round- 
wood harvest volumes. After the logging areas and canopy 
damage intensities were mapped as described in previous 
sections, an equation was applied to convert remotely sensed 
PV change of logged areas to forest canopy gap fraction. The 
equation was derived from intensive measurements of forest 
gap fraction in reduced-impact (low damage) and conven- 
tional (high damage) selective logging areas and co-located 
with Landsat ETM+ satellite imagery processed with the 
AutoMCU algorithm. This general conversion from the 
CLAS-derived PV fractional cover to forest canopy gap per- 
centage was reported by Asner et al. (2005) as: 

Canopy gap=(PV ci ^ s -90.0)/(-0.4) (2) 

PO.ST.pO.Ol 

where PV ciJS and canopy gap are in percentage units. 
P V C l.4s is a planar metric, whereas canopy gap is the hemi- 
spherical canopy opening (S8). The gap -transformed data 
were then used to estimate the volume of roundwood (m 3 ) 
extracted on a per-area basis using an equation drawn from 35 
logging sites in Brazil, Belize, Suriname, Guyana, and Indo- 
nesia (Pereira et al., 2002): 

Wood volume=3.882+108.7(Canopy gap/100) (3) 

P=o.83,p<o.oooi 

Calculated roundwood extraction volumes were then com- 
piled by logging detections (from CLAS), and mean harvest 
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intensities were calculated by dividing the total calculated 
annua] roundwood volume by the harvest area. 
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The foregoing detailed description has been given for 
clearness of understanding only and no unnecessary limita- 
tions should be understood therefrom as modifications will be 
obvious to those skilled in the art. 

While the invention has been described in connection with 
specific embodiments thereof, it will be understood that it is 
capable of further modifications and this application is 
intended to cover any variations, uses, or adaptations of the 
invention following, in general, the principles of the invention 
and including such departures from the present disclosure as 
come within known or customary practice within the art to 
which the invention pertains and as may be applied to the 
essential features hereinbefore set forth and as follows in the 
scope of the appended claims. 

The disclosures of each and every patent, patent applica- 
tion, and publication cited herein including but limited to the 
references listed immediately above are hereby incorporated 
herein by reference in their entirety. This application includes 
and incorporates in its/their entirety the attached CLAS soft- 
ware printouts and supporting documents, all of which are 
individually listed on the ‘Provisional Application For Patent 
Cover Sheet’ accompanying the filing of this application. 

What is claimed: 

1. A method, comprising: 

receiving satellite image data at a computer; 

performing atmospheric correction on the satellite image 
data; 

masking one or more of clouds, water, and deforestation in 
the satellite image data; 

determining fractional material cover at least one of pho- 
tosynthetic vegetation (PV) canopy, non-photosynthetic 
vegetation (NPV) canopy, and bare substrate; 

generating sub-pixel fractional material cover images 
based on the fractional material cover of at least one of 
the PV canopy, the NPV canopy, and the bare substrate; 

performing image differencing to create images of PV 
and/or NPV change; and 

identifying areas of forest disturbance based on the images 
of PV and/or NPV change using pattern recognition. 

2. The method of claim 1, wherein determining fractional 
material cover is performed at a sub-pixel scale from the 
masked satellite image data. 

3. The method of claim 1, wherein the image differencing 
is performed on the sub-pixel fractional material. 


20 

4. The method of claim 1, wherein the pattern recognition 
includes pattern recognition of at least one of the following: 
logging decks, logging roads, skid trails, and tree felling gaps. 

5. The method of claim 1, wherein the step of identifying 
5 further comprises analyzing a cluster of pixels using a moving 

kernel approach. 

6. The method of claim 3, wherein the moving kernel 
approach includes positioning an analysis kernel at a starting 
base kernel with four subset kernels in each of the four cor- 
ners of the base kernel. 

7. The method of claim 6 further comprising the step of 
moving the analysis kernel to its neighboring kernel in at least 
one of the following directions — north, south, east, and west. 

8. The method of claim 7 wherein the step of moving is 
performed if an area of the base kernel tests positive. 

15 9. The method of claim 7 wherein the step of moving is 

performed until all neighboring kernels have been evaluated. 

10. The method of claim 1 further comprising the step of 
performing an audit of probable logging event locations. 

1 1 . The method of claim 10 wherein the audit is performed 
20 to remove potential false logging event locations. 

12. The method of claim 1, wherein the forest disturbance 
includes at least one of the following: agriculture, cattle 
ranching, and a natural hazard. 

1 3 . The method of claim 1 , wherein the step of determining 
25 further fractional material cover comprises spectrally decom- 
posing each pixel of the masked satellite image data. 

14. The method of claim 13, wherein the step of spectrally 
decomposing is performed using a probabilistic spectral sub- 
model. 

30 15. A system, comprising: 

a computer including one or more processors and memory, 
the memory storing instructions for execution by the one 
or more processors, the instructions comprising: 
receiving satellite image data at a computer; 

35 performing atmospheric correction on the satellite 
image data; 

masking one or more of clouds, water, and deforestation 
in the satellite image data; 

determining fractional material cover of at least one of 
40 photosynthetic vegetation (PV) canopy, non-photo- 

synthetic vegetation (NPV) canopy, and bare sub- 
strate; 

generating sub-pixel fractional material cover images 
based on the fractional material cover of at least one of 
45 the PV canopy, the NPV canopy, and the bare sub- 

strate; 

performing image differencing cover images to create 
images of PV and/or NPV change; and 
identifying areas of forest disturbance based on the 
50 images of PV and/or NPV change using pattern rec- 

ognition. 

16. The system of claim 15, wherein the instructions fur- 
ther comprise determining fractional material cover at a sub- 
pixel scale from the masked satellite image data. 

55 17. The system of claim 16, wherein the instructions fur- 

ther comprise image differencing on the sub-pixel fractional 
material. 





